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ABSTRACT 

The cooling theory of neutron stars is corroborated by its comparison with observa- 
tions of thermally emitting isolated neutron stars. An important ingredient for such 
an analysis is the age of the object, which, typically, is obtained from the spin-down 
history. This age is highly uncertain if the object's magnetic field varies appreciably 
over time. Other age estimators, such as supernova remnant ages and kinematic ages, 
only apply to few handful of neutron stars. We conduct a population synthesis study 
of the nearby isolated thermal emitters and obtain their ages statistically from the 
observed luminosity function of these objects. We argue that a more sensitive blind 
scan of the galactic disk with the upcoming space telescopes can help to constrain the 
ages to higher accuracy. 
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1 INTRODUCTION 

The observed thermal states of isolated neutron stars have 
become the primary source to glean useful and interesting 
information about the internal structure of neutron stars 
(NSs). Reconciliation of theoretical cooling curves with ob- 
servations of near by isolated cooling NSs is a challenging 
task (see for e.g. lYakovlev fc Pethickl (|2004l '); iPage et al'l 
l|2006l ) for a comprehensive review). On the theoretical front, 
the problem arises from incomplete knowledge of the compo- 
sition and equation of state (EOS) of matter in the NS core 
at supernuclear densities (p > 10^'*gcm~'^). Many possibili- 
ties can be realized: Depending on the composition the EOS 
can be either soft or stiff, where the stiffness characterizes 
the compressibility of matter, and strongly de pends on the 
intern al degrees of freedom of the system (e.g. Sch aab et al.l 
For example, for a polytropic EOS, P = Kp , a 
larger adiabatic index V yields stiffer EOS. Such an EOS 
generally produces larger maximum masses and radii of NSs 
than its softer counterpart. Softer EOSs can be obtained by 
the introduction of phase transitions in the theory where 
the core of the NS may be composed of boson condensates, 
quark or hyperonic matter. However, the presence of these at 
nuclear densities has been ruled out from t he observation of 
a 1.97±0.04 M0 NS (iDemorest et aLlbOloD . The cooling be- 
havior of NSs depends very sensitively on the composition of 
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the inner core, such that it affects the choice of the neutrino 
cooling process that is do minant in the fir s t 10"^ — 10^ yrs 
of its evolution (see for e.g. lYakovlev et al.l (120011 ) for a de- 
tailed description of all the neutrino emission processes in 
NSs). 

Observations of isolated cooling NSs present its own set 
of challenges in determining their cooling rate. Here, one is 
interested in detecting radiation emanating from the surface 
of the NS which is complicated by the non-thermal emission 
from the magnetosphere. Also, the non-uniform heating of 
the crust due to energetic particles accelerated in the mag- 
netosphere render accurate determination of effe c tive t em- 
peratures hard (see for e.g. the review by lOzell l|2013t ) on 
surface emission from NSs). Also, the unknown composi- 
tion of NS atmospheres leads to the overestimation of their 
temperatures when fitting their spectra with a blackbody 
(e.g. Lloyd et al. (2003)). Apart from the uncertainties in- 
volved in the spectral modelling of NS surface emission, un- 
certainties in their distances can also contribute to poor ef- 
fective temperature and luminosity estimates. Furthermore, 
to place observed isolated sources on cooling curves, accu- 
rate ages are needed that may be hard to obtain as we dis- 
cuss below. 

In this study, we show that, for a large sample, ages 
of isolated thermal emitters can be derived from statistical 
arguments. As the sample size grows, the statistical error 
diminishes. We look at a small group of thermally emitting 
NSs discovered in the ROSAT all-sky survey in Sec. [2] and 
derive their ages statistically in Sec. [3] Lastly, we discuss 
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the possibility of improving the meager sample of isolated 
thermal emitters by detecting more such objects in the up- 
coming eROSITA all-sky survey. 



2 THE MAGNIFICENT SEVEN 

Nearby isolated cooling neutron stars are particularly impor- 
tant for confronting cooling models with observations. They 
were di scovered by Einstein, RQ SAT, and ASCA space tele- 
scopes (iBecker fc Pavlovll2002h . and were further observed 
with the extremely sensitive and high resolution high-energy 
space telescopes Chandra and XMM-Newton. Among all the 
discovered objects, the most interesting are the seven radio- 
quiet thermally emitting isolated NSs (a.k.a. the magnificent 
seven, M7) discovered in the ROSAT all-sky survey (RASS) 
(see for e.g. fe ab crl (2002) for a review). These radio-quiet 
objects radiate predominantly in X-rays with high X-ray to 
optical flux ratios, fx/fopt > 10"*"^. Their soft X-ray spec- 
tra are reasonably well fit by an absorbed blackbody-like 
spectrum with kT < 100 eV and a hydrogen column density 
uh ~ 10^'^ cm^^, indicating small distances d ~ few x 100 
pc. That the thermal emission is coming from majority of 
the stellar surface is confirmed by the small pulse fractions 
< 20% of the X-ray light curves. Spin periods ranging from 

3 - 12 s have been measured for all but one (RX J 185635- 
3754) of the M7 objects (see for e.g. lMereghettil (|201ll '). Table 
[l]). This in conjunction with the measured spin-down rates, 
P ~ 10~^* — 10~^^ s s~^, yields an estimate of the polar 
magnetic field strengths Bp ~ 10^^ G and the character- 
istic spin-down ages Tc ~ 10^ years. Measurements of the 
high proper motions of three of the M7 objects and their 
association thus established to the Sco OB2 complex com- 
prising the Gould Belt yield slightly smaller kinematic ages. 
This discrepancy strongly suggests that the spin-down ages 
are overestimates and the M7 objects in reality are much 
younger. 



2.1 Spin-Down Ages: Poor Age Estimators 

Acco rding to the standard mag netic dipole model of pul- 
sars (jShapiro fc Teukolskvll 19831 ). a rotating NS with a polar 
magnetic field spins down over time by emitting magnetic 
dipole radiation. From the rate of change of the angular 
frequency Q, the spin-down age of the NS can be readily 
determined 
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The spin-down law implicitly assumes that none of the other 
physical characteristics of the pulsar vary over time. This 
may not be the case and, in general, the spin-down law 
l|Lvne fc Graham-Smith. 2006 ) written as the following can 
be allowed to include variation of Bp, the moment of inertia 
7, and the angle between the rotation axis and the magnetic 
dipole axis a, so that 
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where Hi{t) is usually assumed to be a constant and n = 3 is 
the braking index for magnetic dipole braking. Any change 
in n with time naturally yields ages of pulsars that are in 
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Figure 1. Change in the spin-down or characteristic age of an iso- 
lated NS due to various magnetic field decay mechanisms, namely 
the ambipolar (I)rrotational (a = 0.01, q = 5/4) and (S)olenoidal 
(a = 0.15, a = 5/4 ) modes, an d the H all cascade (a = 10, a = 1) 
(See Eq. 2 - 4 of IColpi et ahl 1I2OOOI )). The solid black line de- 
notes t = T, meaning no field decay. This assumes an initial field 
strength Bq = 10^^ G and period Pq = 1 ms. 



conflict with their spin-down ages; Generally, the spin-down 
age should only be taken as a rough estimate to aid in cal- 
culations. 

An independent age estimate is provided by the age 
of the associated supernova remnant (SNR) or massive star 
cluster for younger objects. Establishing such an association 
for older NSs may prove to be difficult since SNRs fade away 
in ~ 60 kyr, and in the same time, due to natal kicks (~ 
500 km s~ ^),NSs may move significantly far away from their 
birth sites l|Frail et al.lll994| ). We plot the spin-down and the 
estimated SNR ages for young pulsars (r < 10^ yrs), central 
compact objects (CCOs), and magnetars (SG Rs and AXPs ) 
along with their timing properties (see for e.g. iBeckerl (|2009l ) 
for a review) in Fig. [1] and it is clear that for most NSs the 
spin-down age is a poor age estimator. The objects that 
have SNR ages smaller than their spin-down ages can be 
explained by having a braking index less than the canomcal 
value, 71 < 3. An excellent example supporting this notion 
is the Vela pulsar which has a very small breaking index 
n = 1.4 ± .2 estimated from an impressive 25- year long 
observation l|Lvne et al.lll99(j ). albeit under the assumption 
that K is still a constant. This yields a spin-down age of 
25.6 kyr, making it appear more than twice as old as its age 
inferred from the standard magnetic braking scenario. This 
result is well supported by the estimated age of t he Vela 
SNR (tsNR - 18 - 31 kyr) l|Aschenbach et al.lll995l ). 

On the other hand, for objects that have spin-down ages 
larger than that of their true ages, that may be inferred 
from their associated SNR ages, it has been argued that the 
magnetic moments decrease in strength over time. There 
are three main mechanisms by which magnetic fields can 
decay in isolated NS, nam ely Ohmic dissipation, ambipo - 
lar diffusion, and Hall drift (jColdreich fc Reiseneggeijll992l ) . 
The timescale over which the fiel d decays substanti aUy due 
to these processes (see for e.g. iHevl fc Kulkarnil ([l99§)) 
determines the dominating process at different stages in the 
evolution of an isolated NS. 

An important consequence of field decay is that it leads 
to an overestimation of the real age of the NS. Following the 
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analysis of (|Colpi et all l|200d ). see Eq. 2, 3, and 4), we plot 
in Fig. [l]the change in the spin-down age of the object over 
time due to the decaying strength of the magnetic field by 
the aforementioned processes. 
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The effect of field decay at late times is apparent from the 
divergence of the characteristic age from the real age. 



3 TRUE AGE ESTIMATES OF ISOLATED 
NEUTRON STARS 

The M7 objects don't have any SNR or massive star cluster 
associations. Therefore, ages for these objects have been de- 
rived from their P and P measurements. In addition, since 
they are nearby objects (d < 500 pc), and due to their large 
proper motions, kinematic ages became a possibility and 
have been estimated for only three of the group members 
(see Table [T]). In this case, one finds that the spin-down 
ages are larger by a factor of 3 — 10 than the kinematic ages. 
Accurate age estimates are extremely important in deter- 
mining the cooling behavior of isolated NSs. Overestimated 
ages used to fit model cooling curves can obscure the deter- 
mination of the true thermal state of these objects. 

3.1 Age Estimates from Population Synthesis 

In the following, we estimate the true ages of the M7 
members by a me thod tha t is motivated by another 
method devised b y Schmidtl (|l968l ) and then applied by 
iHuchra fc Sargend (jl973l ) to calculate the luminosity func- 
tion of field galaxies. The original idea is implemented as 
follows. An apparent magnitude limited sample is first ob- 
tained and it is assumed that the objects in the sample, 
field galaxies for instance, are distributed uniformly in Eu- 
clidean space, such that the luminosity function is indepen- 
dent of the distance. T hen, to each object an accessible vol- 
ume V{M) is assig ned (|Avni fc Bahcal]||l98Gl ). This volume 
depends on the absolute magnitude of the object and gives 
a measure of the volume surveyed for a given object with 
an absolute magnitude M. Basically, V{M) is the maximum 
volume in which the object would be detected had its appar- 
ent magnitude been equal to the limiting magnitude of the 
survey. The whole sample is then divided into bins of size 
dM with objects having absolute magnitude in the range 
[M - dM/2, M + dM/2]. The luminosity function for this 
bin is estimated by adding the inverse of the accessible vol- 
umes for each object in that bin 
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This method provides a non-parametric way of esti- 
mating the luminosity function and it exactly repro- 
duces the true luminosity function within statistical errors 
l|Hartwick fc Schadd[l990l '). Furthermore, this is a very gen- 
eral method which relies on only one underlying assumption 
that the objects are distributed according to their intrinsic 
brightnesses. 

In the case of neutron stars (unlike galaxies) it is rea- 
sonable to assume that the birthrate has been constant over 



the past few million years, and furthermore, the neutron 
stars progress from bright to faint luminosities as the age in 
the same (albeit unknown) way. Under these assumptions 
we can deduce the age of a neutron star of a given 

1 

t{M) ^ - ^{M')dM'dV (5) 
where /3 is the neutron-star birthrate per unit volume. 



3.2 RASS and Population Synthesis 

In the following, we develop a slight variant of the ISchmidtl 
l|l968l) estimator to calculate the true ages of the members of 
the M7 family. The method we develop cannot be completely 
model independent as the distribution of NSs, unlike that of 
galaxies over large scales, is not uniform. Since the progeni- 
tors of NSs mainly reside in the arms of a spiral galaxy, and 
for a natal kick velocity of, say ~ 500 km s~^, the NSs only 
travel a distance of ~ 50 pc from their birth sites within 
~ 10^ years. This is smal l compared to the scale he ight of 
the thin d isk ~ 30 pc jBinnev fc Merrifieldl Il998l ). As a 
result, the ISchmidtl (|l968 ) estimator cannot be used here. 
Instead, we look at the NS progenitor population and cal- 
culate the normalization for each M7 member by counting 
the number of massive OB stars that are found in the acces- 
sible volume Vmax for that object. The population synthe - 
sis method is given in our earlier study jcill fc Hevll(2007l ). 
and essentially requires the assumption of the luminosity 
functio n and spatial distribution of massive OB stars in the 
galaxy (|Bahcall fc Soneira 1980. ). and the distribution of HI, 
which we mod el as a smooth exponentia l disk both radially 
and vertically jFoster fc RoutledgdlioO^ '). 

All M7 objects were discovered by ROSAT which 
scanned the whole sky with a limiting count rate of 
0.015 cts s~^ i n the energy range ~ 0.12 — 2.4 keV (see 
iHiinsch et all (|l999l ) for more details). The complete sur- 
vey covers 9 2% o f the sky for a count rate of 0.1 cts s~^ 
(Vogcs ct all ll999l ) and has yielded the most complete and 
sensitive survey of the X-ray sky. Therefore, it provides a 
perfect flux-limited sample for our study. Next, we calculate 
the weights for each object in the sample by simulating the 
RASS and finding the total number of massive OB stars in 
the volume Vmax, such that the estimated age is given in 
terms of the typical age of their progenitors toB. 
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where Nj is the number of NSs in a small absolute magni- 
tude bin of size dM centered at Mj . Since the sample is of 
marginal size, Nj = 1 in this case. Then, I/Aj.ob gives the 
number of massive OB stars per i^^ object in the sample. 

The ages of NS progenitors are highly uncertain and are 
usually obtained by estimating the main sequence turn-off 
ages of the massive s tar cluster to which the NS may be 
associated (see for e.g. ISmart3 l|2009l ) for a review). Typical 
ages of ~ 3 — 15 Myr hav e been estim ated for the progeni- 
tors of NSs and magnetars. lFiger et al.l (_2005) report the age 
of the cluster of massive stars, containing three Wolf-Rayet 
stars and a post main-sequence OB supergiant, associated 
to th e magnetar SGR 18 06-20 to be roughly 3.0-4.5 Myr. 
Also, iMuno et all (120061 ) report an age of 4 ± 1 Myr for the 
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cluster Westerlund 1 which seems to be the birth site of an- 
other magnetar CXOU J1 64710.2 - 455216. In yet another 
study, IPavies et alj ()2009l ) find the age of the cluster asso- 
ciated to the magnetar SGR 1900 -I- 14 to be 14 ± 1 Myr. 
Since the spin-down ages of magnetars are much smaller 
(~ 10^~* yr) than that of the clusters, the notion that the 
cluster age reflects the age of the progenitor, under the as- 
sumption of coevality of its members, is a valid one. In the 
case of SGR 1806-20 and CXOU J164710.2- 455216, both 
groups find that the progenitor must be a massive star with 
M > 4OM0, except in the last study where the progenitor 
of SGR 1900-1-14 is claimed to be a lower mass star with ini- 
tial MS mass of 17 ± 2Mq. Notwithstanding this last result, 
it has been claimed that magnetars may be the progeny of 
only sufficiently massive stars {PI > 25Mq) (iGaensler et al.l 
I2OO5I ) that would, otherwise, have resulted in the formation 
of a black hole. Although the members of the M7 family 
are endowed with flelds an order of magnitude higher than 
the normal radio PSRs, they are not magnetars and can be 
argued to be the descendants of progenitors not much more 
massive than that of the normal radio PSRs. In that case, 
it is expected that the progenitor age toB will be consider- 
ably longer in comparison to that of magnetar progenitors. 
An upper limit on the ages of M7 progenitors can be placed 
from the age of th e Gould Belt, tos ^ tcs ~ 30 — 60 Myr 
l|Torra et al.l l200dh. 

From Eq. [S]the ages of the sample objects are propor- 
tional to toB, thus significant uncertainty in the progeni- 
tor age will yield erroneous ages. To circumvent this prob- 
lem, we normalize toB such that the estimated age of RX 
J1856.5 — 3754 matches its kinematic age tkin ~ 0.46 ± 0.05 
Myr, which has been measured using its high proper motion 
l|Tetzlaff et aLlbOllI ). This yields a slightly lar ger age for the 
progenitor toB ~ 55 Myr. The corresponding SN rate in the 
Galaxy is ~ 1 per century, in good accord with publishe d 
estima tes, for example 1.9 ± 1.1 reported bv iDiehl et al.l 
l|2006l ). 

Over the last few years, two new candidates have 
been added to the M7 group. Th e first object, IRX S 
J141256.0 + 792204 dubbed Calvera l|Rutledge et al.ll2008l ). 
was actually cataloged in the RASS Bright Source Cata- 
log (|Voges et al.lll999l ) for having a high X-ray to optical 
flux ratio Fx/Fy > 8700. However, its large height above 
the Galactic plane {z ~ 5.1 kpc), requiring a space veloc- 
ity Vz 5100kms~'^, presents a challenge for its interpre- 
tation as an isolated cooling NS like the M7 members (see 
iRutledge et al.l (|2008l ) for a detailed discussion). Also, recent 
X-ray observations of Calvera done with the XMM-Newton 
space telescope found una mbiguous evidence for pulsations 
with period P = 59.2 ms (|Zane et al.l 120111 ). The authors 
of this study argued that Calvera is most probably a CCO 
or a slightly recycled pulsar. The uncertainty in its nature 
fsee lHalpernl (,2011)) doesn't warrant inclusion into our sam- 
ple of radio-quiet is olated NSs. The second object 2XMM 
J104608.7 - 594306 (jPires et aLlbOOgD . discovered serendip- 
itously in an XMM-Newton pointed observation of the Ca- 
rina Nebula hosting the binary system Eta Carinae, appears 
to be a promising candidate (see Table[T]for properties) . This 
object was not detected in the RASS due to its larger dis- 
tance (« 2.3 kpc, based on its association to the Carina neb- 
ula) and higher neutral hydrogen absorption column density 
{Nh = 3.5± 1.1 X 10^^ cm"^). Therefore, the accessible vol- 
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Figure 2. Statistical cooling curve of the nearby thermally emit- 
ting isolated NSs. We plot sample theoretical cooling curves from 
lYakovIev fc Pethidj ||2004| '| for comparison. 



ume Vmax is the ROSAT surveyed volume plus the additional 
volume probed by the XMM-Newton's pointed observation. 

In Table [1] we provide all the relevant data on the sam- 
ple objects including the spectral fit parameters that were 
used to simulate the RASS to obtain Knax. We take the 
calculated ages and plot them against the effective temper- 
atures observed at infinity in Fig. (2] The errorbars on the 
ages correspond to the maximum of the difference in ages 
obtained due to uncertainties in Ttb, Nh (these two param- 
eters are covariant), and the distance. The blackbody tem- 
perature Tbb is obtained by fitting a blackbody spectrum to 
that observed from the source. The temperture, thus, corre- 
sponds to the color temperature of the object and is an over- 
estimation of the effective temperature Te// due to strong 
energy dependence of the free-fre e and bound - free op acities 
of the photosphere (see for e.g. iLlovd et al.l l|2003l )). The 
effective temperature is obtained from Ti,b using a colo r cor- 
rectio n factor fc = Tbb/Teff where 1 < /c < 1-8 (e.g. lOzell 
l|2013 )). For comparison, we als o plot some cooling curves 
from Ivakovlev fc Pethidj l|2004l '). where the non-superfluid 
(No SF) model for a 1.3Mq cannot explain the data. Other 
model curves show NS cooling behavior if proton superfluid- 
ity in the core is taken into account (see lYakovlev fc Pethickl 
(j2004h for more details on the IP and 2P models). 



3.3 Statistical Ages Vs The True Ages 

The ages of isolated NSs have been estimated using different 
methods, namely from the spin-down law, cooling models, 
and kinematics. The method we propose in this study to es- 
timate the true ages of these objects has only been applied, 
in its original form, to estimate the age of white dwarfs from 
their c umulative lumino s ity fu nction in globular cluster (see 
for e.g. iGoldsburv et al.l l|2012f )l. The method itself is purely 
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Table 1. Properties of nearby thermally emitting isolated NSs 



Object 


P 


P 


D 


Tbb 


Nh 


Nob 


T 




(s) 


(10-14 s s-i) 


(pc) 


(eV) 


(lO^o cm-2) 




(Myr) 



RBS 1223^ 
2XMM J104608.7*'' 
RX J1605.3 + 32495 
RBS 1774'' 
RX J0806.4 - 4123" 
RX J0720.4 - 3125^4 
RX J0420.0 - 5022^^ 
RX J1856.5 - 3754' 



10.312 11. 20^ 



20 



9.437* 

11.37^2 

8.39I5 

3.45 
7.06- 



5I8 

r21 



4.1 ± 1.8^ 
5.5 ±3.0^2 
6.98 ±0.0215 
2.8 ±0.31** 
2.97 ±0.07^1 



> 525^ 
2000 

325 - 390^5 
390 - 4301" 
240 ± 25 13 
330+1^° 16 



350 
161 



19 



-18 22 
-14 



118 ± 13 
117 ± 14 
86 - 98 

Q2+19 
^^-15 

78 ± 7 

79 ±4 

t;7+25 
0'_47 

57 ± 1 



0.5 - 2.1 
35 ± 11 
0.6 - 1.5 
4.6 ±0.2 
2.5 ±0.9 

1.3 ±0.3 
1.7 

1.4 ±0.1 



4.5 
0.097 

I. 15 
8.7 
2.9 

II. 5 
0.69 
14.6 



3.839 
2.177 



-fO.958 
0.893 
-t-3.102 
-1.263 

0.631+^9^ 

9 Q/l 7"^" 1 ■'^'^'^ 

U.DO(_Q QQ2 

1.986+J;°f, 

n Q9q-t-0.297 
u.ozo_Q 

7-t-0.059 
-0.045 



0.597"* 



0.014 ±0.003 
0.040 ± 0.026 
0.13 ±0.03 
0.15 ±0.03 
0.23 ± 0.03 
0.26 ±0.03 
0.43 ±0.17 
0.52 ±0.17 
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a statistical one, for which the underlying assumption is that 
the obj ects in the sa mple follow a Poisson distribution (see 
for e.g. llFeltenillQrd )) with a constant production rate. The 
important question to ask here is how good of an estimate of 
the true age is the statistical age. What is the inherent sta- 
tistical error associated to this method of predicting ages? 
Consider the youngest object, RBS 1223, which is also the 
hottest among the eight isolated NSs. The kinematic age of 
RBS 1223 has been found based on its association to pos- 
sible OB associations an d young star clusters to be ~ 0.5 
Myr (|Tetzlaff et al.ll2010 h. although this age should only be 
regarded as the upper limit due to large uncertainties. The 
statistical age of RBS 1223 that we find in our study, given 
that the sample only contains eight such objects, is much 
smaller ~ 0.062^°;°°^ Myr, with uncertainties correspond- 
ing to systematic errors. The statistical error is of course 
much larger than the systematic one. Since the discovery of 
an object in a given volume follows Poisson statistics, the 
relative error scales as l/\fN where A'' is the sample size 
with ages ti ^ t^- Therefore, the error in the age of the first 
object is ±ti where ti is its statistical age. Likewise, the er- 
ror in the age of the eighth object is ±ts / \/8. Evidently, one 
needs a much larger sample to reduce the statistical errors 
to that comparable to the systematic errors. 



4 DISCUSSION 

The ages of isolated NSs are primarily important for con- 
straining their inner structure. In addition, they can be use- 
ful for accurately determining the birthplace of the object, 
if its proper motion is known. Similarly, if the SNR-NS as- 
sociation has been made, then the age of the object can re- 
veal its space velocity and the kinematics of the SNR. The 
knowledge of ages of such objects is also useful for popula- 
tion synthesis models which rely on the spatial and velocity 
distributions, and the birthrates of NSs. In this study, we 
propose a statistical method to estimate the true ages of 
the ROSAT discovered sample comprising the M7. We then 
use the age estimates along with the derived spectral tem- 
peratures to compare the data with some cooling models. 



The strength of this technique, as discussed earlier, lies in 
obtaining a larger sample of coolers. With only eight objects, 
the statistical error is indubitably much larger. 

The statistics can be improved by locating more of these 
objects in th e disk of the G a laxy. I n a recent population syn- 
thesis study, I Posselt et al.l (|2010') find that young isolated 
NSs, that are both hot and bright, with ROSAT count rates 
below 0.1 cts s-i (the ROSAT bright source catalog had a 
limiting count rate of 0.05 cts s'l) should be located in OB 
associations beyond the Gould belt. They also remark on the 
possibility of finding new isolated cooling NSs by conduct- 
ing yet another careful search of such objects in the RASS 
and using the recen tly published XMM-Newton Slew Survey 
(|Esauei et al.ll2006r ). However, they note that ROSAT obser- 
vations are incapable at locating the isolated sources with 
sufficient spatial accuracy, such that many optical counter- 
parts can be found in its large positional error circle. On the 
other hand, although XMM-Newton is much more sensitive, 
albeit with strong inhomogeneities, and can probe deeper 
into the Galactic plane, the slew survey only covers 15% of 
the sky currently. Searching the RASS for new isolated NSs 
may appear to be a promising avenue, however, it is unlikely 
that any new sources will be identified. What is needed at 
the moment is another all-sky survey that is able to surpass 
ROSAT in both sensitivity and positional accuracy. In that 
regar d, the upcoming eROSITA mission (|Cappelluti et all 
l20lj ) shows a lot of promise, and its planned launch in 2014 
makes it very timely. The X-ray instrument eROSITA will 
be part of the Russian Spectrum-Roentgen-Gamma (SRG) 
satellite, equipped with seven Wolter-I telescope modules 
with an advanced version of the XMM-Newton pnCCD cam- 
era at its prime focus. The telescope will operate with an 
energy range of 0.5— 10 keV, a field of view (FOV) of 1.03°, 
an angular resolution of 28" averaged over the FOV, and 
a limiting flux of ~ lO-i* erg cm'^ in the 0.5 — 2 keV en- 
ergy range and ~ 3 x lO^i^erg cm-^ in the 2 — 10 keV 
energy range. The all-sky survey will reach sensitivities that 
are ~ 30 times that of the RASS where the entire sky will 
scanned over a period of four years. 



6 R. Gill and J.S. Heyl 
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